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Abstract 



We present a utilitarian review of the family of matrix groups Sp{2n,'?R.), in 
a form suited to various applications both in optics and quantum mechanics. 
We contrast these groups and their geometry with the much more familiar Eu- 
clidean and unitary geometries. Both the properties of finite group elements 
and of the Lie algebra are studied, and special attention is paid to the so- 
called unitary metaplectic representation of Sp{2n, SR). Global decomposition 
theorems, interesting subgroups and their generators are described. Turning 
to n-mode quantum systems, we define and study their variance matrices in 
general states, the implications of the Heisenberg uncertainty principles, and 
develop a [/(n)-invariant squeezing criterion. The particular properties of 
Wigner distributions and Gaussian pure state wavefunctions under Sp{2n, K) 
action are delineated. 



I. INTRODUCTION 



The symplectic groups form one of the three major families of classical semisimple Lie 
groups, the other two being the real orthogonal family and the complex unitary family [|l|. 
Apart from the groups describing nonrelativistic and relativistic space-time geometries, 
namely the Galilei, Lorentz and Poincare groups, most of the Lie groups encountered in 
physical problems, for example as symmetry groups, belong to either the real orthogonal or 
the unitary families 0. These are multidimensional as well as complex generalisations of 
the rotation group of Euclidean geometry characterising physical three-dimensional space. 
As a result, the intuitive geometrical ideas that go with real orthogonal or complex unitary 
geometries are quite familiar to most physicists. 

It has been realised recently, however, that in several problems both in quantum me- 
chanics and in optics, the real symplectic groups play an important role In the latter 
context, this is so in both classical and quantum theories. More generally, these groups come 
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in very naturally through the canonical formalism of classical dynamics, and its counterpart 
in quantum mechanics. 

Symplectic geometry, on the other hand, differs profoundly from the real or complex 
Euclidean variety [Q. Here the intuitively familiar concepts of length, angle, perpendicularity 
and the Pythagoras Theorem are all absent. In their place we have typically new concepts 
characteristic of canonical mechanics. 

The purpose of this informal and utilitarian survey is to introduce methods based on the 
real symplectic groups to those who are otherwise familiar with the structures of quantum 
mechanics and/or the theory of partial coherence in optics. The intention is to "inform 
rather than astonish", and to describe the main features of the symplectic point of view. 
Our account will not contain complete proofs of all statements presented, but a motivated 
reader should be easily able to supply additional detail and proceed to make practical use 
of these methods. 

The material of this article is organised as follows. Section 2 introduces the group 
Sp{2n, 3?) as the group of linear transformations preserving the classical Poisson Brackets 
as well as the quantum commutation relations among n pairs of canonical variables. In the 
quantum case, the Stone - von Neumann theorem allows us to infer that these transforma- 
tions are unitarily implementable. Section 3 develops some ideas related to real symplectic 
linear vector spaces, specially the concepts of symplectic complement and symplectic rank 
of a subspace, in order to contrast symplectic geometry with real Euclidean and complex 
unitary geometries. Some useful properties of the matrices occurring in the defining rep- 
resentation of Sp{2n,^), and their complex form, are then explained in Section 4. Here 
we also list several useful subgroups of Sp{2n,^), and describe four global decomposition 
theorems - the polar, Euler, pre -Iwasawa and Iwasawa decompositions. Section 5 studies 
the Lie algebra Sp{2n, 9ft), first in the defining representation and then in a general, possibly 
unitary, representation. Convenient ways of breaking up the generators into subsets, and the 
generators of various subgroups, are described. In Section 6 we set up and study the special 
unitary metaplectic representation of Sp{2n, 9ft) and relate it to the generalised Huyghens 
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kernel in any number of dimensions. The characteristic differences between the compact and 
the noncompact generators of Sp{2n, 3fJ) are seen in their dependences on mode annihila- 
tion and creation operators. The former (latter) conserve(do not conserve) the total number 
operator. Section 7 studies the relationship between the metaplectic unitary representation 
of Sp{2n,^), and two often used descriptions of quantum mechanical operators, namely 
the Wigner function representation and the diagonal coherent state representation. While 
the former is covariant, i.e, transforms simply, under the full group Sp{2n, 3?), the latter is 
covariant only under the maximal compact subgroup K{n) = U{n) of Sp{2n,'Si). Section 8 
takes up the questions of defining the noise or variance matrix for any state of an n-mode 
quantum system, both in real and complex forms, and their behaviours under Sp{2n, 3?). In 
Section 9 we carry this analysis further to show that the Heisenberg uncertainty principles 
for any number of modes can be given in explicitly Sp{2n, 3?) covariant forms; a key role 
here is played by Williamson's Theorem relating to normal forms of quadratic Hamiltonians. 
This study leads in Section 10 to the setting up of an U (n)-invariant squeezing criterion for 
n-mode systems. This is the maximal physically reasonable invariance one could ask for in 
these systems, and it can be stated very elegantly in terms of the general variance matrix 
set up in Section 8. Section 11 describes, and motivates, some interesting classes of variance 
matrices with distinctive group theoretic properties, and Section 12 is devoted to a study 
of general centred and normalized pure Gaussian wavefunctions for n-mode systems. The 
transitive action of Sp{2n, 3?) on these wavefunctions, via the symplectic representation, and 
the emergence of a matrix form of the Mobius transformation, are described. Section 13 
contains some concluding remarks. 

II. THE REAL SYMPLECTIC GROUPS Sp{2n, 3ft) 

We consider a classical or quantum canonical system with n degrees of freedom, that is, 
n pairs of mutually conjugate canonical variables. In the classical case these are numerical 
variables written as qr,Pr,r — 1,2, ... ,n. In quantum mechanics we have an irreducible 
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set of hermitian operators qriPr acting on a suitable Hilbert space Ti.. The basic kinematic 
structure is given by Poisson brackets (PB) in one case, by the Heisenberg commutation 
relations (CR) in the other. To express them both compactly and elegantly, we introduce 
the following notation. We assemble the g's and p's into 2n-component column vectors 

i = it) = {qi--- qnPi ■ ■ -Puf, a = 1, 2, . . . , 2n (2.1) 
Then the classical PB's and the quantum CR's are, respectively: 

{^a, ^b} — Pab, 
[L:ib\ = ih/3ab: 

P = iPab) 



^ 1 ^ 

y Inxn Ofjxn J 



(2.2) 



The usual relations stated separately in terms of g's and p's are all contained here; and the 
even-dimensional real antisymmetric matrix /3 will play an important role. 

We assume always that all the ^a, ^-re of Cartesian type: the natural range (spectrum) 
for each of them is the entire real line 3fJ 

Changes of ^, ^ to new quantities ^' given as functions of the old ones such that the 
basic kinematic relations are preserved may be called canonical transformations in both 
situations: 

^' = numerical functions of^ : {C^CU = Pab', 

^' — operator functions of ^ : [^^, = ihPab (2.3) 

Apart from C-number translations (shift of origin) the simplest such transformations are 
the linear homogeneous ones. Each such transformation may be specified by a real 2n- 
dimensional matrix S, the actions being 

S — (Sab) '■ = Sab^b, 

t = S^bib (2.4) 



In either case, the requirements (2.3) lead to a matrix condition on S" [|T|: 

S(3S^ = (3 (2.5) 

This is the defining condition for the real symplectic group in 2n dimensions: 

Sp{2n, ^) = {S = real 2n x 2n matrix | S(3S^ = (3} (2.6) 

The matrix (3 is real, even-dimensional, antisymmetric and nonsingular. It is a "symplectic 
metric matrix". As we see explicitly later, Sp{2n,^) transformations preserve symplectic 
scalar products and the symplectic metric. 

In quantum mechanics the Hilbert space Ti on which the C,a act irreducibly can be de- 
scribed in many ways. The most familiar is the Schrodinger description using wave functions 
on 3?", that is, elements of L^(3ft"). The g,. act multiplicatively while the Pr are differential 
operators: 

n = {^(g) I / d-q\ ^(g) \'< oo}; 

(?r^)(g) = qr^iq), 
d 

(Pr^)(g) = -i'h—i){q), 
dqr - 

g=(gi,...,g„)e3?" (2.7) 

Since the are hermitian and irreducible, and since for any S G Sp{2n, 3?) the transformed 
are also hermitian and irreducible and obey the same CR's, by the Stone-von Neumann 
theorem |^ the change ^ — ^ ^' is unitarily implementable. Thus for each 5* G Sp{2n, 3?) it is 
definitely possible to construct a unitary operator U{S) acting on H such that 

t = Sahib = u{sr'iaU{s), 

U{SyU{S) = 1 on n (2.8) 

This V({S) is arbitrary upto an ^-dependent phase factor. The general composition law that 
follows from the irreducibility of the ia is: 

S', S G Sp{2n, 3ft) : U{S')U{S) = (phase factor dependent on S', S)U{S'S) (2.9) 
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We shall discuss in Section 6 the maximum simplification that can be achieved in this phase 
factor by exploiting the phase freedom in each IA{S). 

III. ASPECTS OF SYMPLECTIC GEOMETRY 

In this Section we develop a few basic concepts related to symplectic vector spaces, so 
that the contrast with Euclidean and unitary geometries can be clearly seen Q. 

Let \^ be a real 2?T,-dimensional vector space, with vectors x,y,- ■ ■. Suppose a nondegen- 
erate bilinear antisymmetric form (.,.) - a "scalar product" - is given on V. Thus for any 
vectors x,y E V, {x,y) is a real number separately linear in x and y; and in addition the 
following hold: 

antisymmetry: {x,y) = —{y,x), 

nondegeneracy: {x,y) = for all y x = (3.1) 

Then is a symplectic vector space, and (.,.) is a symplectic scalar product. 

We have stated the properties of the bilinear form in a basis independent way. It can be 
shown that if in a general basis we express {x, y) in terms of components of x and y in the 
form 

{x,y) = x'^rjy, (3.2) 

involving an antisymmetric nonsingular matrix rj, we can always change to more convenient 
bases in which rj takes on particularly simple canonical, or normal, forms. Two such forms 
are worth mentioning. In one, 7] becomes the matrix j3 of eq.(2.2): 

rj = p : {x,y) = Xiyn+l + X2yn+2 H h a:„7/2„ - Xn+Wl - Xn+2y2 X2nyn (3.3) 

Here the first and {n+l)th components belong to one canonical pair; the second and {n+2)th 
to the second pair; and so on. Another normal form disposes of the canonical pairs one at 
a time : 
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T) — block-diag(i(T2, ia2, ■ • • , i(72) '■ 
{x, y) = xiy2 - X2yi + 2:31/4 - x^ys H (3.4) 

With the normal form (3.3) the meaning of the defining condition (2.5) for symplectic 
matrices becomes geometrically clear: 

SeSp{2n,^),x' = Sx , y' = Sy^{x',y')^{x,y) (3.5) 

In this sense the symplectic scalar product is preserved. Of course this means that we could 
have replaced the condition (2.5) by another entirely equivalent one, using in place of /3 the 
matrix in eq.(3.4). 

Now let us look at linear subspaces Vi C V. In both Euclidean and unitary geometries 
it is well known that all subspaces of the same dimension are basically similar, and cannot 
be distinguished from each other in any intrinsic sense. With symplectic geometry there is 
a difference, as a new concept comes in. 

Given a subspace Vi, we consider the bilinear form {x,y) defined over V, restrict both 
arguments to Vi, and regard the result as a bilinear form on Vi. Now the nondegeneracy 
property may well fail! Thus, there may exist a vector x & Vi such that {x, y) = for all 
y eVi. So we define the rank of this restricted form as the symplectic rank of Vi. 

Symp. rk. Vi — rank {x, y), x and y &Vi (3.6) 

Thus, a Xr, r — 1, ■ • ■ k is a basis for Vi, where k is the dimension of Vi, then Symp. rk. Vi 
is the rank of the k x k antisymmetric matrix {{xr,Xs))- The symplectic rank is necessarily 
an even integer. We have the obvious limits 

< symp. rk. Vi < k (3.7) 

But nondegeneracy of (.,.) over V leads to another nontrivial lower bound, which is effective 
if A; > n: 

(0, 2{k - n))> < symp. rk. Vi < k (3.8) 
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Basically we can say that the symplectic rank of a subspace Vi is twice the number of 
complete canonical pairs contained in V^. Clearly this concept is symplectic invariant. In 
particular two subspaces Vi and V( of the same dimension cannot be mapped on to one 
another by any Sp{2n, 3?) element if they have unequal symplectic ranks. 

As in the Euclidean case, we can pass from Vi to its complement written for convenience 
as V^. But the geometrical significance is quite different. We call the symplectic 
complement of Vi and define it as a subspace of y by: 

V^^^{xeV\ (x, y) = for all y e Vi} (3.9) 

Taking the complement twice gives back Vi. 

(Vi^)^ = Vi (3-10) 

This is as in the Euchdean case. Even the dimensions follow the same rule: 

DimV^^ ^2n-k (3.11) 

This can be shown by using the nondegeneracy of (.,.). But there the similarity ends. It can 
well happen that Vi and V-^^ have nontrivial intersection, and to that extent their sum does 
not give back all oiV. In general, 

V, n v,^ ^ 0, 

Vy^Vi® Vj^ (3.12) 

The two symplectic ranks can be related: 

Symp. rk. V^^ = 2{n - k) + symp. rk. Vi (3.13) 

The extreme case of eq.(3.12), which is very nonintuitive on the basis of Euclidean 
geometric notions, is when symp. rk. Vi vanishes. For this case we give a special name and 
find: 
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Symp. rk. Vi — <^=^ 
Vi is an isotropic subspace of V <^=^ 
(x, y) — for all x,y eVi <^=^ 

ViQV^^ ^k<n (3.14) 

So in this case is contained in V^. The opposite can also happen and then we call Vi a 
co-isotropic subspace: 

Vi is a co-isotropic subspace of V <^=^ 
V-^ is isotropic <^=^ 

\4 D 1/^^ ^ A; > n (3.15) 

So if one of the pair Vi, V^^ is isotropic, the other is co-isotropic. 

An isotropic subspace has dimension k <n, while a co-isotropic one has dimension k >n. 
When they coincide, we have a special situation and name. An n-dimensional subspace Vi 
of V which has vanishing symplectic rank is both isotropic and co-isotropic, and coincides 
with V-^. It is called a Lagrangian subspace. This notion is important in Hamilton- Jacobi 
theory in classical dynamics; it is also relevant in the choice of complete commuting sets of 
operators in quantum mechanics. 

These properties and notions give a feeling for symplectic geometry, and for the ways in 
which it differs from orthogonal and unitary geometries. In particular the notions of length, 
angle and perpendicularity are no longer available. 

IV. PROPERTIES OF Sp{2n, SR) MATRICES, COMPLEX FORM, SUBGROUPS, 

DECOMPOSITIONS 

The matrices in the defining representation of Sp{2n, 3?) obey eq.(2.5). From here many 
useful consequences follow, and we list them: 

(i) Sp{2n, 3fJ) is of dimension n{2n + 1). 
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(ii) P e Sp{2n,^). 

(iii) S e Sp{2n, 3?) ^ -S, S-\ e 5p(2n, 3?), 

(iv) detS = +1. 

(v) 5" e Sp{2n, 3?) =^ eigenvalue spectrum of 5" is invariant under reflection about the real 
axis, and through unit circle (re*^ — > ^e*^); eigenvalues ±1 have even multiplicities. 

(4.1) 

Property (i) can be seen from the number of conditions contained in eq.(2.5), and will be 
confirmed at the Lie algebra level. While properties (ii) and (iii) are easily checked, (iv) is 
rather subtle; an indication will be given later to obtain it. Property (v) is a consequence 
of S and (5*"^)^ being real and related by a similarity transformation. 

Sometimes it is convenient to write S in n x n block form, and then eq.(2.5) becomes a 
set of conditions on the blocks: 



e Sp{2n, 3?) : 

\CD) 

SPS^ = /3 <^ AB'^, CD'^ symmetric, AD'^ - BC^ = l„xn 

S^(3S = /3^ A^C, B^D symmetric , A^D - C^B = l„xn (4.2) 

While it is easy to check (as mentioned above) that S G Sp{2n, 3?) implies G Sp{2n, 3?) 
as well, it is not so easy to pass directly from the first set of conditions above to the second 
set, aside from reconstituting A,B,C,D into S and then passing to 5"^! 

Complex form of Sp{2n, 3?) 

The (3 matrix refiects the precise way in which the real g's and p's have been put together 
in eq. (2.1) to form the 2n component object ^ with real entries. Sometimes it is convenient, 
for instance in dealing with modes of the radiation field, to work with complex combinations 
of the g's and p's - mode annihilation and creation operators defined in this way: 
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% = -^{Qj + iPj),^] = -^{qj - ipj),j = l,---,n (4.3) 
It is useful to arrange these into a new column vector ^^"^^ with non-hermitian entries, 



r(c) 



lai 



1 



n 



V2 



1 -il 



V2 



1 1 

-il il 



(4.4) 



Then the basic commutation relations in (2.2) can be written in two equivalent ways: 



^3 Jab, 



-1 



(4.5) 



Now when we subject ^ to the real transformation 5" e Sp{2n,'^),^^'^^ experiences an equiv- 
alent complex transformation : 



(4.6) 



A + D + i{C - B) A-D + i{B + C) 
^A- D -i(B + C) A + D + i(B -C) J 

Thus 5'^'^) is just a convenient complex form of the real transformation S, much like the 
passage from Cartesian to spherical components of spherical tensors. 

Some Subgroups of Sp{2n, K) 

We shall describe here some useful subgroups of Sp{2n,^). Their dimensions will be 
given, and where it is useful their complex forms exhibited 

(a) GL{n, 9?): This is the n^ -dimensional general real hnear group; in terms of the block 
matrices A,B,C,D it is given thus : 
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(4.7) 



Here the ^'s are subject to a general real linear nonsingular transformation among 
themselves, and then the p's change in a compensating contragredient manner. 

(b) 0(n, K): This is the orthogonal subgroup of GL{n,^), of dimension |n(n — 1). It is 
that part of GL{n, 3fJ) under which the g's and the p's change in the same way: 

A^ D eO{n,^),B ^0 (4.8) 

If we impose the condition detA — +1, we get the subgroup SO{n, 3?) of proper 
orthogonal transformations. 

(c) U{n): Now we come to the n-dimensional unitary group, of dimension n^, a maximal 
compact subgroup within the noncompact Sp{2n, K). We shall sometimes write K(n), 
or simply K, for it. The corresponding symplectic matrices S are identified as follows. 
If we split any U & U{n) into real and imaginary parts we find the properties 

U^X-iY e U(n) , U^U = [/[/t = 1 <^ 
X^X + Y^Y = XX^ + YY^ = 1, 

X'^Y , XY^ symmetric (4.9) 

We can then produce a solution to the matrix condition (4.2)! We find: 

A = D = X^B = -C = Y 



S{X, Y) 



fx y\ 



^-YX, 



e5p(2n,3fJ) (4.10) 



It is an interesting and easy exercise to check the following: If a 2n x 2n real matrix 
is both orthogonal and symplectic, then it is unimodular as well and has to have the 
form S{X, Y) for some U = X — iY & U (n) : 
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0(2n, SR) n Sp{2n, ^) = S0{2n, 5J) n Sp{2n, 5J) 
= K{n) = {S{X, Y)\X -iY e U{n)} 



(4.11) 



The complex form of these matrices is very reveahng: 



'u 0^ 



(4,12) 



So the a's and the a)^s undergo separate unitary rotations, not mixing with one another: 



U e U{n) : a^Ua, ^ U*a) 



(4.13) 



Indeed , K{n) is the maximal subgroup of Sp{2n, 3?) such that a's and a^'s trans- 
form independently. We also have the expected relation between the subgroups 
0{n, 9ft), GL{n, 9ft), and U{n) exhibited above: 



0{n, 9ft) = GL(n,9ft) n U{n) 



(4.14) 



Finally we turn to some Abelian subgroups 0. 

(d) T^: This is a subgroup of dimension ^n{n+l) and may be called the "free propagation" 
subgroup: 



A = D = 1,B = B^,C = 



(4.15) 



The name comes from the actions on q and on p: 



q = q + Bp, p = p 



(4.16) 



Group composition corresponds to adding the B matrices, which explains the Abelian 
nature. 
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(e) T*^'); This is the result of conjugating elements of T^^^ by (3. We call it the "lens" 
subgroup, on account of the action on q's and p's: 

A = D = 1,B = 0,C = C^; 

q' = q, p' = p + Cq (4.17) 

The dimension is again ^n{n + 1), and group composition amounts to adding the C 
matrices. 

Some other subgroups of Sp{2n, 3?) will appear in connection with global decomposition 
theorems. 

Global decomposition Theorems 

Now we describe four useful ways of expressing any S G Sp{2n, 3?) as a product of 
specially chosen factors, either two or three in number. 

(a) Polar Decomposition [0: This says that any S G Sp{2n,^) can be written uniquely 
as the product of two factors, one belonging to the maximal compact subgroup K{n), 
the other to an important sub set n(n) in Sp{2n, 3?). This subset is defined by 

n(n) = {S e Sp{2n, ^)\S^ = S,S positive definite} C Sp{2n, 3?) (4.18) 

and it is definitely not a subgroup. The decomposition reads: 

S G Sp{2n, 3?) : S = S{X, Y)P uniquely, 

S{X,Y)eK{n), P eUin) (4.19) 

Of course by conjugating P with S{X, Y) one could have written the two factors in 
the opposite sequence. The important points here are the global nature of this result, 
and the uniqueness of the factors. From this decomposition one can see that of the 
two possibilities det = ±1 allowed by eqn (2.5) the choice det S = +1 is the only 
one allowed. 
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(b) Euler Decomposition: Next we turn to a decomposition which involves three factors, 
each drawn from a subgroup of Sp{2n, 3fJ), but which is nonunique. Two of the factors 
are from K{n), the third from those elements of 11 (n) which are diagonal and do form 
a subgroup: 

S e Sp{2n, S{Xi, Yi) D{k) 5(^2, Y^), 

5(Xi,yi), 5(X2,r2) G K{n), 

D{k) ^ diag{Ki,- ■ ■ ,Kn,K];;^,- ■ ■ ,K'n^) e n(n), 

Kr>0, r = l,---,n (4.20) 

If one adds the numbers of free parameters in the three factors, one gets the sum 
n{2n + 1) which is just the dimension of Sp{2n,^). Thus the nonuniqueness of this 
decomposition is of a discrete, not a continuous, nature. It stems essentially from the 
freedom to order the first n diagonal elements of D{k) in any way we like. 

(c) Pre-Iwasawa Decomposition: Here we have a three-factor unique decomposition which 
one encounters on the way to establishing the (next) Iwasawa decomposition but the 
factors do not all belong to subgroups of Sp{2n,^). The present decomposition results 
from attempting to reduce the off diagonal block S in a general -S" e Sp{2n, 3fJ) to 
zero, by using an element of K{n) on the right. The result reads: 



S = 



(ab\ 


_{ ^ A 


' Ao 


\ 


ix 








^ A-^ 


I 


.-^ 





X-iY = A^\A-iB), 



Co = {CA^ + DB^)A-^ (4.21) 

Here the matrix Ag is to be chosen symmetric positive definite, and all factors are 
unique. The symmetry of Co A^^ can be checked, so the first factor lies in the lens 
subgroup of eq.(4.17). The middle factor belongs to the intersection CL{n,^) n 
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n(n), which is not a subgroup. And the third factor is from K{n). This particular 
decomposition is of importance in obtaining the generahsed Huyghens kernel, which 
we describe in Section 6. 

(d) Iwasawa Decomposition The polar and pre-Iwasawa decompositions are similar 
in that they involve unique factors, but each factor is not taken from a subgroup. 
The Euler decomposition solves the latter problem, but in the process uniqueness is 
lost. The fourth and last Iwasawa decomposition retains both virtues: it is global, has 
unique factors, and each is taken from a characteristic subgroup of Sp{2n,^). It is a 
result of fundamental group theoretical significance, valid for all simple noncompact 
Lie groups. The three subgroups involved are the maximal compact K{n), a certain 
maximal Abelian subgroup A, and a certain nilpotent subgroup A/". Therefore this 
decomposition is often called the KAM decomposition. We first display it for 5*^(2, K): 

( a 

c d 



S 



ad — be 



1 : 



s 



1 

e 1 



^ ^ Cosip/2 -Sinip/2^ 







'ri/2 



^ = {ac + bd)/{a'^ + b^) 



Sirup/2 Cosip/2 
-oo, oo). 



r] = Ini^a^ + 6^) G (— cxd, oo) 
if = 2arg(a — ib) G (— 27r, 27i] 



(4.22) 



Here the first factor belongs to the subgroup A/", coinciding for n = 1 with the 
lens subgroup T'; the second —r]— factor belongs to the subgroup A; and the third 
—if — factor is from K{1) = 5*0(2). 

For general Sp{2n, 9ft), the situation is more involved. The subgroups A and A/" are: 



A = ^D{k) = diag {ki, ■ ■ ■ , Kn, K^^ , ■ ■ ■ , K^^) | k,. > o| C n( 
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A 

c {A-^Y 



(y ... 0^ 



A 



\ 1 



V 



, A^C symmetric 



J 



> C Sp{2n, ^) (4.23) 



The abelian subgroup A consists of just the elements D{k) that were used in the Euler 
decomposition (4.20). The Iwasawa decomposition for Sp(2n,^) then states that any 
S e Sp{2n, K) can be uniquely expressed as the product of three factors, 



S = 



A 

C {A-Y 



\ 



J 



D{k) S{X, Y) 



(4.24) 



taken respectively from M,A and K{n). The dimensionalities of these subgroups, 
respectively n^, n and n^, add up correctly to n(2n + 1). 



V. THE LIE ALGEBRA OF Sp{2n, 3ft) 

We first study the Lie algebra Sp{2n, 3?) in the defining representation, and then gener- 
alise to any other representation. In keeping with quantum mechanical convention, we shall 
retain a factor of i in the definition, even though this might seem unnecessary in dealing 
with a group of real matrices. 

We examine the form of matrices S e Sp{2n, 3?) close to the identity: 

S = exp {—ieJ) ~ 1 — ieJ, \e\ « 1 : 

SPS^ ^(5^ {(5JY = /3J, {J(5f = J/3, 

r = -J (5.1) 

Thus the generator matrix J is pure imaginary, and both j3J and Jj3 are symmetric. In 
other words in the defining representation we get all possible J's by pre or post multiplying 
all possible pure imaginary symmetric 2n x 2n matrices by (5. Taking the former alternative 
and choosing the simplest possible basis for symmetric 2n x 2n matrices, we obtain the 
following basis for Sp{2n, 3fJ): 
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X. 



(0) 
ab 



XS^\ a,b 



,2n; 



(5.2) 



These matrices can be easily seen to obey the commutation relations 



.(0) 



(0) 



.(0) 



(5.3) 



The structure of Sp{2n, 5ft) is determined by these relations. In a general representation of 

Sp{2n, 5ft) we have generators Xab = X^a obeying 



[Xab, Xcd] — iiPacXbd + f^bcXad + PadXcb + PbdXc 



(5.4) 



Finite dimensional representations of Sp{2n, 5ft) are necessarily nonunitary, hence in them 
the Xab cannot all be hermitian. This is because of the noncompactness of Sp{2n,^). On 
the other hand, in a unitary representation which is necessarily infinite dimensional, we have 

^Ib ^ Xab- 

To help identify the subsets of generators for various subgroups it is useful to use split 
index notation. We use r, s, • • • = 1, • • • , n to label the various canonical pairs; and a, (3, - ■ ■ — 
1, 2 to pick out the q and the p in each pair: 



a,b, 



,2n : a ^ ra, b — > s(3; 

Pab Pra,sl3 ^rs^a/3i 

( , 1^ 



-1 



Then the various components of Xab — Xra,si3 are handled thus: 

Xrl,sl ^rs ^srt 
Xrl,s2 — 



X, 



r2,s2 — ^rs 



(5.5) 



(5.6) 



There are |n(n + 1) V^s, a similar number of Z's, and W^s; in a unitary representation, 
each is hermitian. In this split form the commutation relations (5.4) read: 
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\\Vrg, Vuv\ i^^us^rv ~l~ ^vs^ru) j 

\\Vrsi Zy^y^ — i{5i-uZ sv ~l~ ^rvZsu) 1 
[Vrs, Zuy] = i{6ruWsv + SsuWrv + SryWsu + SgyWru), 

[V,V] =[Z,Z]=0 (5.7) 

Now one can pick out the subsets of generators for various subgroups of Sp{2n,^); we give 
the results in the form of a Table. 

Subgroup Generators 
GL(n, 3?) Wrs 



rs 



SO{n,^) Jrs = Wsr-W, 

U (n) = K{n) Jrs, Qrs = Vrs + Zrs 



rs 



rs 



A Wrr, r — 1, ■ • ■ ,n 

H Wrs for r < s, and all Vrs- 

We have mentioned that any nontrivial finite dimensional representation of Sp{2n, 3?) is 
necessarily nonunitary. It turns out that, with no loss of generality, we may assume that the 
"compact" generators of K{n) are hermitian, while a balance of "non compact" generators 
are antihermitian. That is, in any finite dimensional representation we can assume the 
following: 

Generators of K{n) — compact generators 

— Wrs ~ Wsr, Vrs + -^rs — hermitian; 
Balance of generators = non compact generators 

= Wrs + Wsr, Vrs — Zrs — antihermitian (5.8) 
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The noncompact generators can be arranged into complex combinations with definite tensor 
behaviour under U{n). These combinations are 

Trs Zj.g ^(W^j-s -|- W^^j.), 

Trs = Tsr = Vrs - Z^s + l{Wrs + W^r). (5.9) 

Then the complete set of Sp{2n, 3?) commutation relations (5.4 ) appears in a ?7 (n) adapted 
form P]: 

[^rs;-^iii'] ^ruT sv ^rvT su i 

[T,T] = [T,T] = (5.10) 

We see that Trs and Trs are second rank symmetric tensors under U{n), of contravariant 
and covariant types respectively. While in any representation we can assume Aj^ = Ag^, 
only in unitary representations do we have T^^ = Trs as well. 

VI. THE METAPLECTIC UNITARY REPRESENTATION AND GENERALISED 

HUYGHENS KERNEL 

We saw in Section 2 that for each S G Sp{2n, 3?), on account of the Stone- von Neumann 
Theorem, we can construct a unitary operator ^(5*) such that eq. (2.8) holds. Clearly the 
phase of V({S) is free. We can ask if this ^-dependent phase can be chosen so as to make 
the composition law (2.9) of the Ws as simple as possible. The answer is that this can be 
done, and upon maximum simplification we can achieve 

SuS2e Sp{2n,^) U{Si)U{S2) = ±U{SiS2) (6.1) 

This sign ambiguity cannot be eliminated. So we say that we have here a two-valued unitary 
representation of Sp{2n, 9ft). A more correct or useful statement is that the operators involved 
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provide a faithful unitary representation of the metaplectic group Mp{2n), which is a two- 
fold covering of Sp{2n,^) |T0[. Strictly speaking this means that the argument of U{-) 
should be an element of Mp{2n), not S G Sp{2n, 3?). However, having made this point, we 
shall continue to write U{S) as in eqs. (2.8, 2.9, 6.1). 

The generators of this metaplectic representation of Sp{2n, 3?) are all hermitian; in terms 
of g's and p's they are the quadratic expressions [|ll| 

Vrs = qrQs, Zrs = PrPs (6.2) 

The characteristic differences between the compact and the noncompact combinations be- 
come clear when expressed in terms of a's and a^'s: 

Compact generators : Wrs — Wsr = ^(ajdr — alas), 

Vrs + = {alas + ajor + Srs), (6.3a) 
Non compact generators : Wrs + Wsr = i{alal — dras), 

Vrs ~ Zrs = dial + drds- (6.3b) 

We see that the compact generators of U{n) conserve "total photon number", thus this sub 
group of Sp{2n, 9ft) consists of "passive" transformations. The noncompact generators on 
the other hand do not conserve "photon number", so we may call them "active" generators. 
These properties are expressed thus: 

N = dldr : 

[Wrs-WsrOTVrs + Zrs,N]=0, 

[Wrs + Wsr or Vrs ' Zrs, N] ^ (6.4) 

In fact, single exponentials of i times real linear combinations of the compact generators 
give us operators of the form U{S{X,Y)); while single exponentials of i times real linear 
combinations of the noncompact generators give us operators of the form U{P), P G Tl{n). 
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For this reason the latter may be called "squeezing transformations" |T2| and the polar 
decomposition (4.19) may be read as stating that any metaplectic unitary transformation is 
uniquely the product of a compact passive factor and a non compact active squeeze factor. 
The definition and production of squeezed states are taken up in more detail in Section X. 

The Schrodinger description of the Hilbert space Ti on which the metaplectic represen- 
tation acts has been given in eq. (2.7). In this description, the eigenvectors | g > of the 
commuting position operators appear as a basis: 

qr\q> = qr\q>, 

<q'\q> =5^-\q^-q), 
d 

< q \ Vt = —ifi-T. — < q I (6.5) 
dqr - 

It is useful to know that certain of the operators V({S) have very simple actions on these 
basis vectors. We list them below: 



S{A) 



U ^ 



V 



{A~^f 
U{S{A)) \q) = \ detA 1 1 Aq) 



A e GL{n, \ 



{q I UiSiA)) = I detA l^^/^ ^^-i^ | . 

D(k) = diag^Ki, ■ ■ ■ , k„, k^^, ■ ■ ■ , > : 

U{D{k)) = exp —i ln{Kr)Wrr I , 



(6.6a) 



\ r=l 
/ n \ 1/2 



\r=l 



■ I Kiqi, ■ ■ ■ , f^nqn) 
q I V({D{k)) = (H nr)-'^\n^'qu ■ ■ ■ , ^^n^Qn \ 

r=l 

[-9 l) 



(6.6b) 



^(Ha)) = exp{--grsVrs), 
l^{L{g)) I q) = exp{--q^gq) \ q) 



(6.6c) 
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With the help of these results, and the pre-Iwasawa decomposition for elements of Sp{2n, 3f?) 
described in Section 4, it turns out to be possible to calculate the generalised Huyghen's 
kernel in n-dimensions without too much effort. This kernel is the configuration space matrix 
element < q \ U{S) | g' > of the metaplectic unitary operator U{S). We recall that in the 
case of one degree of freedom, this kernel has the following form [|rT| : 



s 



{q I UiS) I q') 



I 



c d 

g-i7r/4 



e Sp{2, 3?) 
exp [i{dq^ - 2qq' + aq'^)/2hb] , & 7^ 0; 



1 1/2 



(6.7) 







I exp{icq^ /2a)5{^ - q')/ 
( These results are, strictly speaking, valid only if S is sufficiently close to the identity, the 
point being that W(-) is actually a representation of Mp{2) and so carries as argument an 
element of this group. Therefore the generalised Huyghens kernel is not expressible totally 
in terms of 5* G 5*^(2, 3ft)). We may regard the case 6 7^ as generic. This generalises nicely 



to any number of dimensions and we find: 



S 



{q I UiS) I q') 



A B 
C D 

-jn7r/4 



G Sp{2n, 



detB ^ : 



exp 



I 

2h 



(6.8) 



detB I 

{q^DB-\ - 2q'^B'\ + q'^B'^Aq'} 

The nongeneric case when detB = has to be handled carefully - then the kernel collapses 
to a lower dimensional expression with a certain number of delta function factors. However, 
having given an indication of the structure involved, we will not go into any further details. 



VII. Sp{2n, K) ACTIONS ON WIGNER AND DIAGONAL COHERENT STATE 

REPRESENTATIONS 

Let r be any quantum mechanical operator, specified in the Schrodinger representation 
by its configuration space kernel {q \ T \ q'). We can see that if T is conjugated by U{S) for 
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general 5* G Sp{2n, 3fJ), the change in the kernel involves an integral transformation in which 
the generalised Huygens kernel (6.8) and its complex conjugate both appear. We can ask 
whether there is any other way of specifying or describing f such that this change takes a 
simpler form, not requiring any integrations at all. Indeed there is, and it is given by the use 
of the techniques due to Weyl, Wigner and Moyal (WWM) We describe this aspect, 
and then go on to another practically important way of describing operators, namely via 
the diagonal coherent state representation, and its behaviour under Sp{2n,^). 

Let us hereafter set ^ = 1. From the configuration space kernel (q' | f | q') of f we obtain 
its Wigner distribution or WWM representative by a partial Fourier transformation : 

f ^ W{0 = (27r)-" I rq'iq - ig' | f | g + p')exp{tq' ■ p), 
{q\t\q') = J dyW{^{q + q'),p)exp{-tp.{q - q')) (7.1) 

Here W{C,) is a function on the classical phase space corresponding to the quantum system, 
with arguments which are 2n classical c-number g's and p's. As seen above, one can recover 
the operator F from its WWM representative W{^) unambiguously. Then one finds that 
under conjugation by the metaplectic operators U{S), the changes in F are very simply 
expressed in terms of W{^) [|T^: 

f' = U{S)-^tU{S) ^ W'iO = W{SO (7.2) 

This behaviour of VF(0 f^ict be regarded as the key or characteristic property of 

the WWM method in quantum mechanics; we may say that this description of operators is 
covariant under the full symplectic group Sp{2n, 9ft). 



Next we turn to the diagonal coherent state description of operators F [|I6|. For n degrees 
of freedom the coherent states are defined as usual by 

f 1 " "1 
\ z) =exp< --Y.\ P +J2^rdl > I 0), 

I ^ r=l r=l J 

dr\0) = 0, 

dr \ z) = Zr \ z), Zr & C,Z = (zi, • • • , Z„) (7.3) 

25 



These are normalised states, no two being orthogonal, and can be written also as the result 
of phase space displacement operators acting on the ground state | 0): 



n 



:) = exp < ^{zrdl — z*dr) > | 0), 



.r=l 



{z' I z) = exp j-^-^'^-^' - ^^^z + ^'^^} (7-4) 
The coherent states form an (over) complete set; the resolution of the identity 

i = /n^u)(2i (7.5) 



shows that any vector | tp) can certainly be expanded using them 

" d^Zr 
r=l 

ijj{z) = {z* I ip) = exp{—-z^z) (entire analytic function of z) (7.6) 

Moreover the overcompleteness allows expansion of any operator T in the form of an integral 
over projections to these states : 

f= ff[^<Piz)\z){z_\ (7.7) 

For given F, this expansion and the weight function (j){z) are unique, however the latter 
could in general be a distribution of quite a singular kind. 

Now we look at the behaviour under Sp{2n,^). It turns out that the states | z) have a 
simple behaviour only under the "passive" maximal compact subgroup K{n) C Sp{2n,^): 

U = X-iY e U{n) : 

U{S{X,Y)) \ z) =\Uz) (7.8) 

This can be traced to the fact that the generators of K{n) involve only terms of the form 
a^a, as we see in eq.(6.3 a). On the other hand, the effect of U{P), for any P G Il{n), on 
I z) involves an integration using a suitable kernel, namely the generalised Huyghens kernel 
expressed in the coherent state language [0, which for the Sp{2, 3?) case can be written in 
a simple form in terms of complex SU{1, 1) parameters A fi: 
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IC{z,z';S^) = {z'\U{S)\z) 



1 



1 + 



J" 



X 



^2|/m(^-A)| 
exp {(A* + z^ + z z'* + (A'^ - z'*^} 
Im {fj, — X) , i.e. ,ip 2mi ; 



2|Re(A + /i) 
1 + (A + /i) Re (A + /i) """"P 



Re (A + /i) + z^*)^ 
1 + (A + /i) Re (A + /i) 



Im (/i — A) = , i.e. ,ip = 2mT , n = —1, 0, 1, or 2 . 



S 



A /i 



\ = ^{a + d + ic — ib) 
fi = ^{a — d + ib + ic) 



(7.9) 



(Here each individual square root is defined to have a positive real part). 

In contrast to the WWM result (7.2), we now have covariance under K{n) alone: 

r = u{s{x, Y)y^t u{s{x, y)) ^ ^'(z) = 0(t/s) (7.10) 

Active elements of Sp{2n, 3?) change in a manner involving a nontrivial integral transfor- 
mation. 

VIII. QUANTUM NOISE MATRICES AND THEIR 5p(2n, 5ft) BEHAVIOUR 

Let p be the density operator of any (pure or mixed) quantum state. For simplicity alone 
let us assume that the means of vanish: 



(U) = Tr{pU) = 



(8.1) 



Non zero values for these means can always be reinstated by a suitable phase space displace- 
ment. The variance or noise or second order moment matrix of the state p is then defined 



as follows 12 : 
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V = (Kfe) = 

Va, = V,a = \{{LL}) 



(^2)„ = ^{{qr,Ps}), 

(Vsls = iPrPs) (8.2) 

This is a real symmetric 2nx2n positive definite matrix subject to further matrix inequalities 
which express the uncertainty principles (see below). 

We must note the compact way in which we are able to express Vab as a phase space 
integral involving the WWM representative W{C,) of the density operator p. This is a 
consequence of the general rule [0] 



Tr(pexp {iCPi}) = I d'-^ W{0 exp {tCP^} (8.3) 

valid for any numerical ^o- Indeed, one can regard this property of W{C,) as being as basic as 
the symplectic transformation rule (7.2). A particular case of eq.(8.3), when = 0, shows 
that W{C,) is normalised, because Trp = 1: its phase space integral is unity. But we must 
remember that though for hermitian p the function W{C,) is real, it may not be nonnegative 
everywhere. 

The change in the noise matrix V when p is changed by a symplectic transformation is 
now easily obtained by exploiting eq.(7.2) along with the above expression for Vab in terms 
of W{^). We have the extremely simple law, for any S G Sp{2n, 3?): 

p' = U{S)pU{S)-^ = SVS'^ (8.4) 

(Note that in comparison to the change f ^ f ' in eq.(7.2), S has been replaced by 
S^^ here). We may say that V undergoes a symmetric symplectic transformation, which 
preserves its symmetry and positive definiteness. 
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The information contained in V can also be given in complex form using second order 
moments of a's and a^'s: 



yB* A* ^ 

yic) _ yic)* _ l/icic) tWtU. 



■rs — ■^sr 



ar, a 



Bf-g Bgy ijXj-Qjg 



(8.5) 



Thus is a hermitian 2n x 2n positive definite matrix subject to further uncertainty 
inequalities given below. 

The relations connecting these two forms of the noise matrix are: 

Vi = ]^{A + A* + B + B*}, 
V2='-{A-A*-B + B*}, 

V^^]-{A + A* -B-B*} (8.6) 

Here the matrix VL is as given in eq.(4.4). And in place of eq. (8.4) we have the equivalent 
transformation law 



(8.7) 



IX. WILLIAMSON'S THEOREM AND UNCERTAINTY PRINCIPLES 

For general 5" e Sp{2n, 3?) the transformation law (8.4) for the noise matrix is not a 
similarity transformation, it is a similarity only if S* = S{X,Y) G K{n). Normally we 
expect that the diagonalisation of V will require a matrix belonging to the group S0{2n, 3fJ). 
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However, the fundamental Theorem of Wilhamson comes to our rescue |T8| . This theorem 



is a complete answer to the question: Given a real symmetric 2n x 2n matrix V , what is the 
maximum simplification we can achieve in the form of V' = SVS^ by allowing 5* to vary all 
over Sp{2n, 3fJ)? For general V, the normal or canonical form of V is not a diagonal form; 
however the Theorem shows that in case V is positive (or negative) definite to begin with, 
then we can certainly choose an S so that V is diagonal. By suitable further rescaling and 
ordering of elements we can then achieve the following: 

V real symmetric positive definite: 

SVS'^ = diag{Ki, ■ ■ ■ , ki, ■ ■ ■ , k„), 

/^i < /«2 < ■ • • < f^n, 
Suitable 5 e 5p(2n,3fJ). (9.1) 

We shall call this the Williamson normal form of V . In general, the are not the eigenval- 
ues of V at all. Also note that when V is in this form, then we have V^^^ = V. 

In the Williamson normal form we see that for each canonical pair qr,Pr we have equal 
uncertainties: Aq^ = Ap^ = k^^. Also all the off-diagonal variances vanish. Therefore the 
complete statement of the uncertainty principles for all degrees of freedom would be 

Kr> 1/2, r = l,---,n (9.2) 

A given 2n x 2n real symmetric positive definite matrix V is quantum mechanically realisable 
as the noise matrix of some state p if and only if in its Williamson normal form (9.1) every 
diagonal entry is greater than or equal to one-half. This is an Sp{2n, 3?) invariant statement, 
to express these uncertainty principles directly in terms of V without actually passing to its 
Williamson normal form [|I^]. There are several ways of doing this and for illustration we 
quote just one, expressed both in terms of V and V^^^: 

V + -(3 = hermitian positive semidefinite , 
y*-'^-* -|- -E3 = hermitian positive semidefinite. (9.3) 
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We emphasize that any V, V^'^^ obeying these conditions are quantum mechanically realisable 
- they are necessary and sufficient. 

There is a subtle distinction between the matrix V being brought to diagonal form, and 
the matrix V^^^ being diagonal. From the relations (8.6) among the two sets of submatrices 
we can see easily: 

V diagonal -^^^ V\ , V3 diagonal , V2 = , 
<^=^ B real diagonal 

^ l/^'Miagonal; (9.4a) 
V^^^ diagonal <^=^ A real diagonal , B = 
=^ V\ , VJj diagonal, V2 = 

V diagonal (9.4b) 

This means that the set of states which have diagonal V'^'^^ is a subset of the set of states 
having diagonal V: the former is a more restrictive condition, so fewer states obey it. In 
the particular situation when V is in Williamson normal form, it is not only diagonal but 
in addition Vi = V3; thus, in spite of the general statement (9.4 a), we do find V''''^ also 
diagonal in that case. In complex form, the relevant consequence of Williamson's Theorem 
for V'^'^^ is that with the transformation law (8.7), V^'^^ can be " diagonalised" . 



X. A U(N)-IN VARIANT MULTIMODE SQUEEZING CRITERION 



The noise matrix Y for any state of a single mode system is two dimensional and has 
the form 



V 



V 



10.1) 



A(g,p) (Ap)2 

with obvious meanings for the various matrix elements. The usual Heisenberg uncertainty 



relation (with h = 1) which reads |T9| 
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AqAp>^ (10.2) 



can, as is well known, be strengthened to the statement [BO 



deW ^ {Aqf{Apf - {A{q,p)f > 1 (10.3) 

This is in fact the content of the conditions (9.3) in this case, and this is Sp{2, 3?) invariant. 

The state p with variance matrix V is usually said to be a squeezed state if either one 
of the two diagonal elements of V (but of course not both) is less than one half. However 
this definition possesses no useful or interesting continuous invariance at all. A definition 
of squeezing possessing invariance under the maximal compact U{1) or 5*0(2) subgroup of 
5*^(2, 3?) is this: the state p is squeezed if and only if the lesser of the two eigen values of 
V is less than one half. A state squeezed in the former sense is squeezed also in the latter 
sense but not conversely, so the former is more restrictive: 

Aq or Ap < —j= =^ lesser eigenvalue of ^ < - (10-4) 
V 2 2 

The U{1) - invariant squeezing criterion has been used in several studies. It is of course 

clear that we cannot ask for any more invariance in the squeezing criterion, for example it 

would be meaningless to think of an 5*^(2, 9ft) invariant squeezing criterion. 

Motivated by the above, we now describe a criterion for squeezing for states of n mode 

systems [jl2| . Suppose that for a given state p the noise matrix V already has some diagonal 



element less than one half. Then we say that the state is manifestly squeezed. However it 
may happen that every diagonal element Vaa of V exceeds or equals one half, yet squeezing is 
buried and not manifest. We have divided elements of Sp{2n, 9ft) into passive U{n) elements 
and active H(n) elements. We would like to have a definition of squeezing invariant under 
passive U (n) transformations. Such a definition is the following: 

p is a squeezed state <^=^ {S{X, Y)VS{X, Y)^)^^ < - , 

some X — iV & U {n) , 
some a = 1, 2, ■ ■ ■ , 2n. (10.5) 
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Thus, p is squeezed if V shows it to be manifestly so, or if this happens after a suitable 
passive transformation. 

We know that in general V cannot be diagonalised by similarity transformations within 
K{n) - we expect to have to use matrices from the much larger group S0{2n, 3?), which may 
be noncanonical. Nevertheless it is remarkable that the above U{n) - invariant squeezing 
criterion can be expressed in terms of the eigenvalue spectrum of V: 

p is a squeezed state <^=^ 

1{V) = minimum eigenvalue of \^ < - (10.6) 

This is precisely equivalent to the definition (10.5) 

Obviously the squeezed or nonsqueezed nature of a state p is unchanged by passive ele- 
ments of Sp{2n,^) lying within K{n), since then V undergoes a similarity transformation 
which leaves /(V") unaltered. To change /(V^),and so the status of a state, in either direc- 
tion, we must use active noncompact elements lying in n(n), if we wish to do so within 
the framework of Sp{2n,^) transformations. This justifies our using the term "squeezing 
transformation" for the metaplectic unitary operator V({V) for elements P G n(n). 



XI. SOME INTERESTING FAMILIES OF VARIANCE MATRICES 

Let us denote by S the set of all allowed noise matrices V, i.e., all physically realisable 
ones obeying the uncertainty principles: 

S = = 2n X 2n real symmetric positive definite \ V + positive semidefinitej (11.1) 

This is an n{2n + 1) parameter family. We have seen that a general G iS is not 
diagonalisable by elements within K{n). We ask: can we characterise the subset Sk C S 
consisting of all noise matrices which are diagonalisable using elements of K{n) ? The 
answer is that this can be done rather elegantly, and for this the complex form V^^^ of V is 
more convenient |T^: 



33 



<Sk = {V E S\AB = symmetric} C S : 
V eSk^ S{X, Y)VS{X, Yf = diagonal, 
someU = X -iV eU{n) (11.2) 

The subset Sk is an n(n + 2) parameter family. 

There are two further subsets of Sk which are interesting [|T2|. We call them the "her- 
mitian" family Sh, and the "Gaussian" family So - they are of dimensions and n{n + 1) 
respectively. The definition of Sh is motivated by the fact that the general transformation 
rule (8.7) for V^^^ becomes very simple if we have a K(n) element: 

UeUin): V(^) — ^ S^'\U)V^'^S^'\uy, 
A — y UAU-\ 
B — > UBU^. (11.2) (11.3) 

So if i3 = to begin with, it remains zero; and A being hermitian can be diagonalised by 
some U . This would then result in both y^'^) and V becoming diagonal. Thus, including the 
uncertainty conditions (9.3), the definition of Sh is : 

Sh = \v e Sk \ B = Q , A-]^-l positive semidefinitej d Sk ^ S (11.4) 

One can verify that no member of this family Sh is squeezed. 

Next to the family Sg- Here the definition involves the "noncompact" subset n(n) of 
Sp{2n,^) : 

Sg = {V eSK\2V e n(n) C Sp{2n, ^)} d Sk d S (11.5) 

It is a fact that such noise matrices are physically realisable, i.e., they obey the uncer- 
tainty conditions, and are diagonalisable within K{n). Further, except for V" = | ■ 1, 
every other V E Sq is squeezed. So we see incidentally that this is the only common ele- 
ment in Sh and Sg- All in all, we have: 

^H, Sg C Sk C S 

5Hn5G = |i-i} (11.6) 
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XII. GAUSSIAN PURE STATES, GAUSSIAN WIGNER DISTRIBUTIONS 



The generators of the metaplectic unitary representation of Sp{2n, 3?) are quadratics in 
g's and p's. The most general centred n-mode Gaussian wave function involves a quadratic 
in the g's in the exponent. It turns out that the former act in very nice and compact ways 



on the latter We describe the main features briefly in this Section. 

A general centred and normalised n-mode Gaussian pure state can be parametrised by 
two real symmetric n x n matrices u and v, of which the former is positive definite: 



:i2.i) 



For {u,v) = (1,0) we get the ground state of the isotropic oscillator in n-dimensions: 



V'(i,o)(g) = TT "^^exp(-^g^g), 

dr ^(1,0) = 0. 



[12.2) 



The calculation of the WWM representative for il'{u,v){q), and of the noise matrix, are easy 
since only Gaussian integrals and moments are involved. The results are: 

- r 1 1 

W(^u,v){0 = (2vr) "y d''q'ij(^u,v){q - p')^{u,v){q + -q'Texp{iq' ■ p) 



7r-"exp{-^^G(M,t;)^} 



G{u, v) 



u + vu vu 



/ 

V 

{S{u,v)-YS{u,v) 



-1 



U V u 



S{u, v) 



V{u,v) 



1 



1 



e Sp{2n, 3?) ; 
S{u, v)S{u, f )^ G Sg 



:i2.3) 



While the calculations leading to these results are elementary, it is worth paying attention to 
the structures involved. The WWM representative W(^u,v){0 is expected to be a Gaussian, 
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with a positive definite parameter matrix G{u, v) in the exponent. What is interesting is 
the factorization of this matrix in terms of an Sp{2n, 3fJ) - matrix S{u, v); an added feature 
is that this S{u, v) is an example of the product of the first two factors in the general pre- 
Iwasawa decomposition (4.21) for any S e Sp{2n, 3?)!. That the noise matrix V{u, v) should 
be essentially the inverse of G{u,v) is clear from eq.(8.2); it is then the product structure 
for G{u, v) that results in V{u, v) being an element of the family Sq defined in the previous 
Section. Incidentally we also see that n(n) C Sp{2n,^) defined in eq.(4.18) and used in 
defining the family Sq in eq. (11.4) can be described more explicitly using S{u,v): 

U{n) = {S{u,v)S{u,vf\S{u,v) e Sp{2n,^) 

u and V real symmetric n x n matrices, u positive definite} (12.4) 

We can argue further along similar lines, concerning the form to be expected for the action 
of an operator U{S) on any 4'(u,v){q)- From eq. (7.2) it is clear that W(^u,v){0 must get 
mapped onto another Gaussian; the fact that it arises from a pure state wave function must 
also be retained; and the transformation rule (8.4) plus the explicit form of V{u, v) in (12.3) 
means that any iP{u,v){q) can be mapped onto '4'{i,o){q) by a suitable U{S). All this is indeed 
true. We find that 

'ipM = (phase factor )W(5'(M,t'))V^(i,o) , 
V(u, v) = S(u, v)V{l, 0)S(u, vf, 

V(l,0) = i-1 (12.5) 

So Sp{2n, 3fJ) acts transitively on the set of Gaussian pure states. Since the subgroup of 
Sp{2n, 3?) leaving the particular state ^^(1,0) invariant (apart from phases) is just K{n) - 
dr annihilates "0(1,0)) a-nd the generators of K{n) in the metaplectic representation are of 
the form aja^. More precisely, the situation may be described as follows. One can easily 
establish that the behavior of '^(1,0) under K{n) is given by: 

U = X-iY e U{n) : 
U {S{X, Y)) ^^(1,0) = VDetU V'a.o) (12.6) 
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The sign ambiguity here explains the appearance of the metaplectic group Mp{2n). The 
stabihty group of ip{ifl) - the subgroup of Sp{2n, 9ft) leaving this vector strictly invariant- 
is thus the n"^ — 1 parameter subgroup SU{n) G U{n). Correspondingly the orbit of ^/'{i,o) 
under Sp{2n, 3ft), made up of the vectors W(S')'?/'(i,o) for all S G Sp{2n, 3ft), is an {n{n + 1) + 
1) parameter family. It consists of the vectors e^'^ip(u,v) for < a < 27r and all allowed 
(«,!>). This orbit is essentially the coset space Sp{2n,^)/ SU{n); there is a one to one 
correspondence between vectors e''°'4'{u,v) and points in this space. Suppressing the phase 
we can next say that the set of density matrices '4'(u,v)'4'lu,v) ; equally well the set of 
representative vectors ip(u,v) , is essentially the coset space Sp{2n,^)/U{n). 

Turning next to the effect of ^(5") on ip(u,v), we see that apart from a phase factor it has 
to result in ipiu'y) for suitable u' and v'. The formula for this change is a beautiful one: 

U{S)i)i^u,v) = (phase factor) 

A = (m - w)"^ ^ A' = {iu' - v')~^ 

= {A A +B){C A +D)-^ (12.7) 

The Gaussian WWM function in (12.3) arose from a pure state. Suppose now we consider 
a general centred normalized Gaussian phase space distribution with a general parameter 
matrix G 

^^^(0 = vr-"(detG)^/2exp(-^^G0, 

G = real symmetric 2n x 2n positive definite matrix. (12.8) 

The question is: when is this a WWM function corresponding to some pure or mixed 
quantum state? Thanks to Williamson's Theorem, the answer is elementary [|15|. The noise 
matrix would clearly be 

V = ^G-' (12.9) 

So we start from the given Gaussian 14^^(0 > pass to its V, then go to the Williamson normal 
form (9.1) ofV : 

37 



V diag(Ki, • • • , Kn, «i, • • • , «n), (12.10) 

and then demand that each Kr be greater than or equal to one half. This is a complete 
necessary and sufficient condition for WcH) to be a bonafide WWM phase space distribution; 
as we have seen, however, this condition can be stated directly without going to the normal 
form: 

Wg{C) is a WWM distribution <^=^ + ip — hermitian positive semidefinite (12.11) 

XIII. CONCLUDING REMARKS 

In this review we have tried to convey the main features of the family of real symplectic 
groups Sp{2n, 3?), and have outlined some problems in optics and quantum mechanics where 
they are useful. Our account has been descriptive and suggestive, omitting detailed prop- 
erties of various statements made. We believe that any interested reader wishing to apply 
symplectic techniques to any concrete problem would be well equipped for the purpose, and 
able to supply necessary details. 

Some general remarks - partly to counter apparently common misconceptions - may 
be useful at this stage. The group Sp{2n, K) comes in when we define linear canonical 
transformations on given canonical variables. It is the fact that the commutation relations 
and hermiticity are maintained that is responsible for the existence of unitary operators U{S) 
implementing these transformations. In particular, U{S) is unitary whether S e K{n), when 
S itself is unitary, or S ^ K{n), for example S e II(n) in which case S is hermitian rather 
than unitary. Correspondingly, the 'Hamiltonians' generating these unitary operators U{S) 
are hermitian quadratics in the g's and p's. For the case of the most general unitary evolution 
via a general hermitian Hamiltonian, the action on ^'s and p's (equivalently on a's and a^'s 
) definitely preserves the canonical commutation relations and hermiticity properties, but 
the transformed operators may not be linear combinations of the original ones. 
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The maximal compact subgroup U{n) of Sp{2n, JJ) has naturally played an important 
role in our considerations. The ra-mode squeezing criterion described in Section 10 has a 
built-in U{n) invariance. As a result, for a state with a given variance matrix V, squeezing 
if present may be manifest (one of the diagonal elements of V is less than 1/2) or may be 
hidden (this happens only after a suitable U{n) transformation). This makes it clear that 
our squeezing criterion is weaker than the usual one stated directly in terms of the diagonal 
elements of V, since these would rule out the hidden case. Correspondingly there are more 
states which are squeezed by our criterion than by the usual one, which in any case has 
much less invariance built in. 

The metaplectic group is essential to describe properly the stucture of the unitary op- 
erators implementing linear canonical transformations on the canonical variables (in this 
respect the notation U{S) with S G Sp{2n, 3?) is inadequate - the argument in U{S) should 
be an element of Mp{2n)). The importance of this group is seen in, for example, the cal- 
culation of geometric phases for cyclic evolution of squeezed states, the interpretation of 



the Guoy phase, etc ||2^. It has also been shown elsewhere that the metaplectic group is 
relevant in setting up operator Mobius transformations for one degree of freedom |T7| . 

Work on extending U{n) invariant notions beyond quadrature squeezing is in progress 
and will be reported elsewhere. Thus, for two-mode systems, one can classify squeezing 
transformations into U{2) invariant equivalence classes, study bunching and antibunching 
with such invariance etc [1^. These notions can be generalized to n mode systems as well. 
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